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Abstract: We study the Green's function of a gauge invariant fermionic operator in a 
strongly coupled field theory at nonzero temperature and density using a dual gravity de- 
scription. The gravity model contains a charged black hole in four dimensional anti-de Sitter 
space and probe charged fermions. In particular, we consider the effects of the spin of these 
probe fermions on the properties of the Green's function. There exists a spin-orbit coupling 
between the spin of an electron and the electric field of a Reissner-Nordstrom black hole. On 
the field theory side, this coupling leads to a Rashba like dispersion relation. We also study 
the effects of spin on the damping term in the dispersion relation by considering how the spin 
affects the placement of the fermionic quasinormal modes in the complex frequency plane in 
a WKB limit. An appendix contains some exact solutions of the Dirac equation in terms of 
Heun polynomials. 
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1 Introduction 

Through gauge/gravity duality [1-3], a charged spinor field in an asymptotically anti-de Sit- 
ter (AdS) spacetime in a classical limit can be used to model strongly interacting fermions 
in field theory. While the start of this program can be traced back to refs. [4, 5] which solve 
the Dirac equation in pure AdS spacetime, with refs. [6-8] there has been a resurgence of 
interest in the subject focused on fermions at nonzero charge density in the hope of modeling 
strongly interacting cousins of Fermi liquids. These so-called non- Fermi liquids are believed 
to underly some of the interesting physics of heavy fermion compounds and high temperature 
superconductors. Initial gauge/gravity duality studies focused on charged black hole back- 
grounds. By tuning the parameters of the fermionic field, both Fermi liquid and non-Fermi 
liquid behavior can be obtained [9]. 

These holographic models of strongly interacting fermions appear to be delicate to con- 
struct. The charged black hole is subject to a wide variety of potential instabilities in the zero 
temperature limit. For example, if the field theory contains an operator dual to a charged 
scalar field in the bulk, the black hole can develop scalar hair in a holographic superfiuid 
phase transition [10-13]. If a four fermion interaction term with the right sign is added to 
the Dirac Lagrangian, there can be a Bardeen-Cooper-Schrieffer phase transition at low tem- 
peratures in the bulk [14]. Even with no extra terms in the Lagrangian, if the fermions have 
a large enough charge, the condensation of a Fermi sea in the bulk will modify the geometry, 
producing an electron star at low temperatures [15-17]. The field theories dual to these elec- 
tron stars exhibit the usual Fermi liquid behavior. In contrast, the non- Fermi liquid behavior 
found by [9] occurs in the limit where the Fermi sea outside the black hole is very small. 
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The delicate nature of these systems aside, they seem to be promising toy models to 
address some of the questions surrounding strongly correlated electron systems. The current 
paradigm surrounding these toy models appeals to ideas of confinement in large N gauge 
theories [18-21]. The charge of the black hole should be carried by deconfined degrees of 
freedom ~ non-gauge invariant fermions behind the horizon. The added Dirac field is dual to 
a confined degree of freedom, i.e. a gauge invariant fermion or mesino. Just as in QCD where 
the mesons and baryons interact weakly with each other in a large N limit, these mesinos 
form a Fermi liquid which by definition is effectively weakly interacting. It is the conjectured 
fermions behind the horizon which lead to non- Fermi liquid like behavior. 

Our goals in this paper are modest. We would like to provide a more careful considera- 
tion of spin physics and spin-orbit coupling in these holographic systems than has appeared 
heretofore in the literature. A qualitative discussion of spin-orbit effects appears in ref. [22] in 
the context of coupling fermions to a d-wave holographic superconductor, but we shall try to 
be more quantitative and precise here. By spin orbit coupling, we mean that a bulk charged 
fermion moving perpendicular to an applied electric field experiences an effective magnetic 
field that splits the degeneracy between the two spin states. 

Let us start with an electrically charged black hole in AdS4^ at nonzero temperature to 
which we add a spinor field. Using the standard gauge/gravity duality dictionary, we may 
compute the quasinormal mode (QNM) spectrum of the spinor field which will allow us to 
deduce where the retarded Greens function for the dual gauge invariant fermionic operator 
has poles [23, 24]. When the charge of the black hole is large enough, the imaginary parts 
of many of these quasinormal modes will be small, and we plot in figure la the real part of 
the freqency of these modes versus momentum. One may think of these curves as dispersion 
relations for fermionic quasiparticles in the field theory. Alternately, one may think of these 
curves as locations where the field theory may have a nonzero density of states. (To compute 
the actual density of states, we need the full Green's function including the residues of the 
poles, and these residues may vanish at special points.) The system is rotationally symmetric, 
and one can envision the full k dependence by rotating the graph around the w-axis. The 
details of this numerical computation are presented in section 3. 

The similarity between figure la and figure lb is one of the central observations of this 
paper. Figure la resembles four copies of figure lb. Figure lb is the dispersion relation for 
a nonrelativistic two dimensional electron gas with a spin orbit (or Rashba) coupling. The 
Rashba Hamiltonian can be written 



where A is the Rashba coupling constant, fi a chemical potential, a the Pauli matrices, nies 
the electron effective mass, and z the unit vector perpendicular to the gas. 

From the bulk spacetime point of view, the similarity between these two figures is straight- 
forwardly explained. Electrons with the dispersion relation (1.1) can be produced by a two 
dimensional slab-like geometry with a strong electric field perpendicular to the slab. On the 



H = 



2mcs 



Xa ■ {z X k) — fj, , 
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gravity side, the charge of the black hole provides the electric field. The slab lies between 
the boundary of AdS on one side and the horizon on the other. More precisely, deriving a 
Schrodinger equation for the spinor field, one finds a potential barrier between the well in 
which the spinors live and the horizon. Tunneling through the barrier produces the small 
negative imaginary part of the QNMs. 

For the 2+1 dimensional field theory, this similarity naively presents a puzzle. The usual 
derivation of the Hamiltonian (1.1) is intrinsically 3+1 dimensional and relies on the presence 
of the electric field and a notion of spin-orbit coupling. In the context of heavy fermion 
compounds and strange metals, one anticipates that spin should be essentially an internal 
SU{2) symmetry of the electrons. There are no strong magnetic or electric fields in these 
compounds, and the Fermi surface or surfaces should be spin degenerate. There is no obvious 
mechanism for breaking the SU{2) symmetry of these strongly interacting fermions at nonzero 
density. 

The solution to this puzzle is that by the rules of the AdS/CFT correspondence, the 
dual field theory is intrinsically relativistic. Spin is not an internal symmetry but instead 
implies a corresponding transformation rule under the Lorentz group. In 2+1 dimensions, 
angular momentum dualizes to a scalar. Massive free fermions satisfying the Dirac equation 
(7^p^ — im)^ = carry a spin determined by the sign of their mass ^ sgn(m) (see for example 
[25]).-^ In contrast, massless fermions, because the little group is too small, carry no spin at 
all. We believe that the fermions in our field theory are massless for two reasons. The first is 
that the field theory is conformal and a mass term breaks scale invariance. The second is that 
while a mass term for the fermions breaks parity in 2+1 dimensions, the state we consider 
in the field theory appears to be parity invariant. In support of this claim, note that the 
Hamiltonian (1.1) is invariant under the parity operation that sends {kx,ky) — )• {—kx,ky) and 

There is an alternate intrinsically 2+1 dimensional way of motivating the Hamiltonian 
(1.1). Our black hole background is dual to a conformal theory at nonzero chemical potential 
fj, and temperature. The presence of energy and charge density identifies a preferred Lorentz 
frame = (1,0,0). From an effective field theory point of view, it is natural to expect a 
modified Dirac equation of the form [26] 



where p^^ = {uj,kx,ky) is the four momentum and F and G are arbitrary functions of uj and 



k = yjk"^ + ky. Note we have not included a bare mass in this expression. Choosing gamma 
matrices 7* = ia^, 7^ = ax and 7^ = ci^ and setting F and G to zero, we recover (1.1) without 
the k"^ term and with A = 1. (In other words, the Rashba coupling itself is the Hamiltonian 
for a massless relativistic fermion in 2+1 dimensions.) To add the k'^ term, we may posit that 
G{oj, k) ~ k"^ which is allowed by the symmetries. To get the different bands in figure la, 
we may additionally posit the existence of several species of massless fermions with different 

^In our conventions, the gamma matrices obey {7*^, 7"} ~ 2ri^'^ where 77^" — { — h +). 



[(1 + F)p^ + + G)u,,] 7^^^ = 



(1.2) 
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Figure 1. (a) The dispersion relation in the boundary theory for a bulk fermion with charge times 
chemical potential divided by temperature iiq/AirT — 25 and mass times AdS radius mL = 2. (b) 
The dispersion relation for a fermion with a Rashba type coupling. The Fermi surface in both cases 
is indicated by the dashed line. 

charges qi, replacing fi with fiqi in (1.1). Ideally, we would like to derive this effective Dirac 
equation from an action, but we do not know how. 

One may object on technical grounds to the nonextremal black hole background used to 
produce figure la. From earlier work on the electron stars [15-17], for the large charge, low 
mass fermion chosen, it is clear that the bulk Fermi sea will cause a back reaction of the black 
hole background. In principle, one should use the numerically computed electron star metrics 
to produce figure la. However, the phase transition between the charged black hole and the 
electron star is third order [16], and there are many qualitatively similar features between the 
charged black hole and electron star backgrounds. Using the numerical electron star metrics, 
which are only exact in an Oppenheimer-Volkoff approximation, should not change the story 
in a qualitative way.^ 

We begin this paper by revisiting the Dirac equation in these charged black hole and 
electron star backgrounds in section 2. We show, using the Pauli-Lubanski pseudovector, 
how to identify the different spin components of the fermion. Next in section 2.1, we review 
how to compute the QNMs of a spinor field in these backgrounds, and connect the QNMs 
to poles of the retarded Green's function in the dual field theory. Section 2.2 reviews how 
to convert the Dirac equation into an effective Schrodinger equation for the spinor for use in 
a WKB approximation. The WKB limit will give us qualitative insight into the nature of 
this spin-orbit coupling. Finally, in section 3, we employ the machinery set up in the earlier 
sections to compute the QNMs of the spinor in a charged black hole in AdS 4, giving the 
details behind figure la. Section 3.1 contains a discussion of the numerical solution of the 
Dirac equation, while 3.2 discusses the Dirac equation in a WKB limit. The WKB calculation 
contains some unusual features which we discuss at length. We are able to show how the WKB 

■^See [19, 27] for progress going beyond the Oppenheimer-Volkoff approximation. Note that in the 
Oppenheimer-Volkoff approximation, the electrons are heavy and the spin splitting is very small. 
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approximation gets the sign of the imaginary part of the quasinormal modes correct. Making 
use of Heun polynomials, appendix A contains some exact analytic solutions of the Dirac 
equation for a charged spinor in a black hole in AdS^. 

Note Added: After this paper was completed, we became aware of [28] which has overlap 
with ours. 

2 The Dirac Equation Revisited 

To have a gauge invariant fermionic operator O^r in a dual field theory, we consider a spinor 
^' in a curved spacetime with the action 

Sq, = -i j d'^+^x^/^'^ij^'Df, - m)^ , (2.1) 

where \I' is a spinor of mass m and charge q, and ^' = ^'^^7-. Vielbein indices are underlined, 
and related to coordinate indices by 7- = e-^7'^. The covariant derivative is 

D^ = df, + ^a;^,a67- - iqA,, , (2.2) 

where uj is the spin connection, 7— = ^[7-, 7-], and is a gauge field. 

With AdS/CFT applications in mind, we make the following simplifying assumptions on 
the metric and gauge field. We assume a translationally and rotationally invariant metric of 
the form 

ds"^ = gttdt^ + gxx{dx^ + dy^ + . . .) + Qzzdz^ , (2.3) 

where the diagonal metric components g^f_i{z) depend only on a radial coordinate z. Addi- 
tionally, we assume that At is the only nonzero component of the vector potential and that 
it is a function only of the radial coordinate z. In other words, there is a radial electric field 
whose strength may vary as a function of the radial direction. The type of spacetimes we 
have in mind are Reissner-Nordstrom (RN) black holes and electron stars in AdS, both of 
which obey this set of assumptions. 

Given these assumptions, the Dirac equation takes a particularly simple form. The 
spin-connection term in the Dirac equation can be eliminated by using the rescaled spinor 
'tjj = {—gg^^)^^'^'^. The equation of motion for tp is 

[7^(a^-zgA^)-m]^ = 0. (2.4) 

Translational symmetry in the directions orthogonal to z suggest taking a Fourier transform 
of V': 

i'ixf^,z) = I ^e*'=''-'>fc(z), (2.5) 

where = (w,k). Because of rotational symmetry, we assume without loss of generality 
that the spatial momentum is in the x direction. By plugging a single Fourier mode ip ~ 
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-iuit+ikx 



^l^{z) into eq. (2.4), we obtain the equation of motion for i/j: 



i\/-5**7-(w + qAt) + \/g^j-dz + i^fcf^-i-k - m 



V' = 0. 



Specializing now to d = 3, we make the following choice of gamma matrices 



7- 



io"2 
ia2^ 



7- 



(T3 

. as, 



7- 



ai 
-fJi 



—iai 
iai , 



(2.6) 



(2.7) 



This choice is consistent with the conventions in [29] and allows for an easy generalization 
to a spacetime of arbitrary dimension. Because we have set momentum in the y direction 
to zero, the Dirac equation for the four component spinor ij: decouples into equations for 
two-component spinors ip = {ip+,ilj-)'^: 



\/-g^^a2{uj + qAt) + ^fg^azdz ± i^f(f^a\k 



m 



V'i = 0. 



(2.8) 



We argue that iIj± correspond to fermions with opposite spin. To see the spin direction, 
consider the Pauli-Lubanski pseudovector 



Wa 



1 



^abcd J P ) 



(2.9) 



where J— = ^[7-, 7-] and P- = —iD-. We wih show that ij) = (^+,0) and = (0, V'-) are 



eigenstates of Wy. Acting on e 



-iojt+ikx 



ip{z), P- = —ie-^gf^^Di, is given by 



P- = (\/^(w + qAt), y/^k, 0) . 



(2.10) 



Thus we obtain W., 



/g^^k + \J —g^^{w + qAt) i\fg^dz 





\ 



iVg^dz ^g^^k + \J —g^*-(w + qAt) 

fg^k — iy + qAt) —i\fg^dz ) 



(2.11) 



We find then that the Dirac equation (2.8) can be written in terms of the Wy component of 
the Pauli-Lubanski pseudovector: 



1 



1 
-1 



(2.12) 



which means that the spin of V'+ is in the y-direction, and the spin of V'- is in the opposite 
direction. Because we assumed that the momentum is in the x-direction, and the system 
has rotational invariance, the direction of spin is always perpendicular to the plane defined 
by the momentum and the radial direction z. The Dirac equation decouples into these spin 
eigenstates. 

This decoupling of spin eigenstates is in good agreement with flat space, nonrelativistic 
intuition [22]. Starting with a massive fermion moving in the j;-direction, the fermion in its 
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own rest frame will experience both an electric field in the radial z-direction and a magnetic 
field in the y-direction. This magnetic field will induce an energy splitting between the 
y-spin up fermion and the y-spin down fermion. The WKB limit we consider later will 
further strengthen this non-relativistic intuition. There is also a possible coupling between 
the curvature and the spin. 

2.1 Prom quasinormal modes in the bulk to dispersion relations in the boundary 

In this section, we would like to take advantage of the well known AdS/CFT relation between 
QNMs of the bulk spacetime and poles in the Greens functions of the dual field theory [23]. 
To that end, we begin by discussing the QNM boundary conditions for the Dirac equation 
(2.8). 

We will assume that our metric in the limit z — )• is asymptotically of anti-de Sitter 
form: 

1 

gtt ~ gxx ^ gzz ^ ^ ■ (2-13) 

There are two linearly independent solutions of the two component Dirac equation which 
approach the boundary as 




XI ~ Q J and X2 ~ I . I . (2.14) 

We apply the "Dirichlet" boundary condition that the X2 solution vanish. Using these bound- 
ary conditions, the relation between the scaling dimension of the dual operator and the 
mass of the spinor is A = m + d/2 [4]. The unitarity bound on the scaling dimension restricts 
us to m > —1/2.'^ 

As the "Dirichlet" boundary condition is Hermitian, the quasinormalness of the modes 
comes from the boundary condition applied in the interior of the geometry. We assume that 
gtt —7- at some Zh > 0, and at this horizon, the phase velocity of the spinor wave function is 
in the positive z-direction. For example, for a non-extremal blackhole, we may take the time 
and radial metric components to vanish as —gtt ~ g^^ ~ 4:TrT{zh — z) where T is the Hawking 
temperature. In this case, our ingoing boundary condition is 

^±^{zr,-zr^l^-^(^^. (2.15) 

One may also consider more exotic situations, for example the Lifshitz geometry in the interior 
of the electron star [15]. 

Solving the Dirac equation with "Dirichlet" boundary conditions at z = and ingoing 
boundary conditions at the horizon is only possible for a discrete set of complex frequencies a; 
called QNMs. As we tune k in the Dirac equation, a given QNM will trace out a curve in the 



^Using the "Neumann" boundary condition, we may consider spinors with dimension A — — m + d/2 and 
< 1/2. 
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complex uj plane. Provided the imaginary part of the QNM is small, this curve is essentially 
a dispersion relation for a quasi-stable particle. 

To relate the QNMs to poles in the fermionic two-point function in the dual field theory, 
let us briefly recall how to compute these two-point functions. The first step is to solve 
the Dirac equation with ingoing boundary conditions at the horizon and arbitrary boundary 
conditions at z = 0, yielding a solution of the form -i/^-t = h±xi+a,±X2 for each spin component 
where a± is interpreted as proportional to a source in the dual field theory and h± as an 
expectation value. The retarded Green's function is in reality a 2 x 2 matrix, but our choice 
of momentum has diagonalized it. Through the theory of linear response, the retarded Green's 
function is = ida/sba/afs [29]. If vanishes, Gr will have a pole. 

2.2 A Schrodinger form for the Dirac equation 

While at least three papers have considered the WKB limit of the Dirac equation in this 
AdS/CFT context [9, 17, 18], these papers have largely ignored the effects of the spin of the 
electron. Our focus here shall be on the spin. 

The WKB approximation in this case assumes that the dimensionful parameters m, u, 
qAt, and k are large compared to the scale dzln^± over which the wave function varies. 
We may capture this limit by introducing a small parameter h multiplying dz in the Dirac 
equation. 



and then expanding in h. 

Our first step is to convert the Dirac equation into Schrodinger form. As noted in this 
AdS/CFT context by [9], there is no unique procedure. Given an equation in Schrodinger 
form, one has the usual freedom to reparametrize the coordinate y = f{z) and rescale the 
wave function (p — t- Zcj) subject to the constraint Z'^f is a constant. However, there is an 
additional functional degree of freedom associated with converting a Dirac equation into 
Schrodinger form; one can introduce the scalar wave function </)-i- such that 



The condition that (f)± satisfies a Schrodinger type equation —d1(f)± + V(f)± = puts three 
constraints on the functions a, /3, 7, and 5, but leaves a family of potentials V{z) parametrized 
by an undetermined function [30]. 

One simple choice is to select (f)± to be proportional to the first component of the two 
component spinor ip±: 




(2.16) 




(2.17) 




(2.18) 



where we have introduced a normalization factor 




(2.19) 
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Given these substitutions, (p± satisfies a Schrodinger equation of the form 



- 9>± + V±k,m(l>± = , (2.20) 

where the potential function is 

V,,m{z) = ^ {g^.m' - Z^,Z,) - ""f^^^^^^'^ + ^^'^ • (2-21) 

Note that the spin dependence of the potential enters only at subleading order. At leading 
order in ^ is independent of the sign of k. 

At leading order in ^, the spinor potential V^.m is exactly what one obtains for a charged 
scalar particle in this curved spacetime. If we start with the scalar wave equation [D^D^ — 
m^)^ = and let $ = e-*'^*+*'''=Z(z)</>(z) where Z = ./<h~z{-g)~^''^ , then we obtain -dlcf) + 
Vgcj) = Q where 

Vs{z) = p {g,,n? - Z_kZk) + Zdl^ , (2.22) 

and we have introduced factors of h analogously to the spinor case. 

There is an important difference between the scalar and the spinor case. The normaliza- 
tion factor Zk may vanish at a point in spacetime where the energy of the particle is equal 
to the local chemical potential plus a k dependent correction. At such a point, the spinor 
potential Vk^m will have singularities at subleading order in h. We will see in the next section 
how these subleading singularities mean that while scalar QNMs can lie in the upper half 
plane, the spinor QNMs will not. 

Before getting into the details, we can say something about the number of zeroes Zk 
possesses in the interval < z < z^. In general, we will associate the boundary value of 
At with a chemical potential: A.t{0) = fi; at the horizon we set At{zh) = 0. Given the 
boundary conditions described in section 2.1 for the metric g^u, close to the horizon we find 
that Zk ~ ijj/ATTT{zh — z). At the boundary, we find instead that Zk{0) = u + qfi — k. Thus, 
if uj + qfi — k and lv are of the same sign, there will be an even number of zeroes between the 
boundary and the horizon. If they have opposite sign, there will be an odd number. In the 
cases we consider below, > k — to, and so the parity of the number of zeroes is determined 
by the sign of cj, odd if u < and even if cj > 0. 

Before continuing to an example, we would like to point out another simple choice of 
Schrodinger equation. We can let (p± be proportional to the second component of -i/'i rather 
than the first. To compute this alternate Schrodinger equation, we do not need to do the 
calculation again. Note instead that eq. (2.16) is invariant under k — )• —k, ni — )• — m, 
and ip± — )• o^l^p±. Thus the Schrodinger equation for the second component is given by 
replacing Vk^m with V-k-m- Interestingly, this symmetry implies that the potentials Vk^m and 
V-k-m have the same quasinormal mode spectrum (being careful to exchange the boundary 
conditions on the two components of 'ip± as well). The existence of this second isospectral 
potential allows for some cross checks when we perform a WKB analysis below. 
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3 An Example: AdS-Reissner-Nordstrom Black Hole 



To study a fermionic system at nonzero temperature and density, we use the RN black hole as 
the background geometry, solve the Dirac equation coupled to a U(l) gauge field, and obtain 
the QNMs. We will solve the Dirac equation both numerically and using WKB. Our main 
interest is the location of the poles of the Green's function. Most of our results are for the 
m = 2 spinor in AdS4^. 

We consider a charged black hole in AdS4 in the Poincare patch for which the metric has 
the form 



ds' = ^ \^-f{z)dt' + dx' + dy^ + j^j . (3.1) 

The conformal boundary is at z = 0, and the horizon is at z = Zfi, where f{zh) = 0. We set 
L = 1 and z^ = 1, which are allowed by two scaling symmetries [13]. The charged black hole 
solution and its probe limit are as follows: 

2 2 

RN-^d54: ■^ = ^" (l + x)'^^^^^^ ^^^^ At = K'^-z); (3.2) 

probe limit: f = 1 — z^ with At = /u(l — z) . (3.3) 

We work in the probe limit fj, = and fiq = held fixed. In this limit, the electric field 
does not backreact on the geometry. If we use the full solution to the RN black hole, the 
qualitative features will not change, as we will discuss at the end of section 3.1. 

3.1 Numerics 

To solve the Dirac equation (2.8) numerically in this spacetime, we approximate the ingoing 
solution by a Taylor series near the horizon z = 1, integrate numerically to a point near the 
boundary z = 0, and fit the numerical solution to the boundary expansion ■0-1- = b±xi+CL±X2- 
The numerical integration was performed using Mathematica's NDSolve routine [31]. As we 
are interested in QNMs and the corresponding Green's function singularities, we focus on the 
value of the source a± . 

Figure 2 presents a density plot of |o-|-| as a function of complex uj. We see that the poles 
of the Green's function in the complex w plane have the following features: 

• As we increase the chemical potential /Xg, there are more and more poles along the 
negative real uj axis. Decreasing fig moves these poles to the right. Once they cross the 
imaginary axis, they begin to move away from the positive real axis. See figure 2. 

• If we fix |Uq and increase the momentum k, the poles move to the right. As a; = 
corresponds to the Fermi surface, when a pole crosses the imaginary axis, we obtain a 
Fermi momentum kp. See figure 3. 

(n) 

We can obtain many Fermi momenta kp and many Fermi surfaces. The poles close to the 
real lo axis correspond to quasibound states, and the number of them at = equals the 
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;i,=5,k=0 A',=15,k=0 ^^=25,k=0 




-10 -5 5 lU -10 -5 5 10 -10 -5 5 10 



Re(ijj) Re(w) Re(w) 

Figure 2. Motion of tlie poles in the complex uj plane as we increase fig. From left to right, fig = 5, 
15, 25. Other parameters are m = 2 and fc = 0. 
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Figure 3. Motion of the poles in the complex u plane as we increase fc. From left to right, fc = 5, 
10, 15. Other parameters are m = 2 and /i^ = 25. 

number of Fermi surfaces. The number of Fermi surfaces grows linearly with the chemical 
potential, which can be seen by comparing figure 2 and figure 11 in the next section (or 
equivalently tables (3.20) and (3.21)). 

For the other spin (k — )• —k), the poles will be in slightly different locations. The spin 
splitting can be seen more obviously by looking at figure 1. The figure is a plot of the 
dispersion relation and was constructed by superposing density plots of |a+| and |o_| in the 
k-io plane. The density plot of ja+l corresponds to the parabolas on the left and |a_| to the 
parabolas on the right. 

For comparison, we also plot the dispersion relation for a massless m = bulk fermion in 
figure 4. (Here only |a+| is plotted. The density plot for the other spin component |a_| is the 
mirror image.) Unlike the massive case, there do not seem to be well defined quasiparticles 
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Figure 4. Dispersion relation for a massless fermion with Hq = 20. Its mirror image k — > —k is for 
the opposite hehcity. There are no poles in the Green's function at fc = 0. 
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Figure 5. The left plot shows the Fermi surfaces with spin splitting. The parameters are m — 2 and 
fiq = 25. The right plot illustrates the directions of spin for massive bulk fermions. 

close to the A; = axis. For the m = case there are no poles in the Green's function at = 
and m = 0. This fact can be checked analytically by solving the Dirac equation to obtain 
Gr{uj, k = 0) = i [9]. From a bulk point of view, the absence of these poles is presumably 
related to the absence of a rest frame for a massless particle. As can be seen in figure 4, 
poles do appear for k ^ 0. For the field theory dual, the interpretation is more obscure. Ref. 
[27] relates the absence to strong interactions with a background continuum of states existing 
inside an "IR lightcone". Why then the interactions are suppressed for larger bulk masses 
still needs to be explored. 

Figure 5 (left) shows the location of the Fermi surfaces in momentum space and also 
demonstrates the spin splitting effect. As we discussed above, the spin of the bulk spacetime 
fermions is perpendicular to the momentum and the electric field. Thus the bulk Fermi 
surfaces are spin polarized, as shown schematically in figure 5 (right). A similar effect has 
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T/p=Q.l T/^=0.Q1 T/Ai=Q.0Q01 
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Figure 6. Temperature dependence of the dispersion relation with fixed /i^ = 25. This is the only 
plot that uses the full solution of the RN black hole in this work. The chemical potential for the 
three plots are /z « 1.7665, 3.2218, 3.4615 (from left to right). The extremal case corresponds to 
/I = 2-\/3 ~ 3.4641. These are spin up fermions. 



been studied in spin Hall systems [32] and observed in experiments [33] for electrons that 
while confined to a plane still have a three dimensional spin. 

Note that for the 2+1 dimensional field theory fermions, the spins shown in figure 5 (right) 
are misleading. From a purely 2+1 dimensional perspective, we argued in the introduction 
that the fermions are both massless and spinless. The extra degrees of freedom producing the 
second Fermi surface come from the hole states. The dispersion relation for the hole states 
has been bent upward, emptying out the infinite Fermi sea and producing a second Fermi 
surface. 

Before moving on to WKB, we promised a discussion of the validity of our probe approx- 
imation. The temperature for the charged black hole solution eq. (3.2) is 

If we lower the temperature by increasing the chemical potential n, more Fermi surfaces will 
appear, but the size of the outer Fermi surface will remain roughly the same, as figure 6 shows. 
By comparing figure 6 and figure la, we can see that if we consider the high temperature 
regime, i.e., T//i > 0.1, there is no essential difference between the full RN solution and the 
probe limit in eq. (3.3). 

3.2 WKB 

The rough picture of the WKB analysis of the Schrodinger equation (2.20) with the spinor 
potential function Vk^m (2.21) is easily explained. Schematically, we may write our potential 

Vk,m = J^ + ^ + V2. (3.5) 
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Figure 7. The leading order potential Vq for the spinor or scalar when uj ^ 0. There exist quasi- 
bound states when Wmin < W < Wmax- 

Considering only the leading order term Vq = gzz^n^ — Zi^Z_i^, there is a barrier at the 
conformal boundary z = provided m? > 0. At the horizon z = 1, the potential Vq is 
unbounded below provided u)'^ > 0. For intermediate values < z < 1 and appropriate 
choices of k, and q, one may find a potential well where Vq < separated from the horizon 
by a barrier where Vq > (see figure 7). For discrete choices of to, the wave function will 
satisfy a Bohr-Sommerfeld type quantization condition, and the potential well will support 
quasi-bound states. Tunneling through the barrier to the horizon gives w a small imaginary 
part. Note in the special case a; = shown in figure 8, the potential at the horizon becomes 
a barrier, suggesting that the quasinormal modes have very small imaginary part close to the 
origin. 

The detailed picture of this WKB analysis is more intricate. At leading order in h, 
the Schrodinger potentials for the spinor and scalar are identical, and the results will be 
independent of the spin that is this paper's main focus. Moreover, there is an important 
qualitative difference in the QNM spectrum for the spinor and scalar that is not captured 
at this leading order. The charged scalar QNMs may lie in the upper half of the complex u 
plane signaling a perturbative instability [11, 12], while the spinor QNMs will not. Yet, as 
we will see, the leading order WKB analysis would put the spinor QNMs in the upper half 
plane as well. 

To capture the effects of spin, we will keep the first subleading term in the h expansion 
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tj=0, k>kF 

0)^0, k<kF Vo(z) 




Figure 8. The leading order potential Vq for the spinor or scalar when uj — 0. There exists a bound 
state provided k is not too large. 



of Vk^m- Our WKB wavefunction is then 



<^WKB = exp (±i IVVo[l + dz^ . (3.6) 

The classical turning points are defined in terms of the zeroes of Vq. 

To set up the WKB problem, let zi < Z2 be the turning points bounding the classically 
allowed region. Similarly, let Z2 < z^ bound the potential barrier. We begin by writing formal 
expressions for the WKB wave functions to the left and right of the three points zi, Z2, and 
Z3 valid to next to leading order in h: 



^0 



) 



Figure 9 portrays a typical potential with the regions labeled in which the six WKB wave 
functions are valid. 

The wave functions in adjacent regions will be related by connection matrices Mj such 
that tij+i = MiVi where Vi = {Ai,Bi)'^. Given these connection matrices, we can obtain a 
semi-classical quantization condition on uj by applying boundary conditions. We will take 
Dirichlet boundary conditions at z = that Ai = 0. At the horizon, we have ingoing 
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Figure 9. A typical spinor or scalar potential at leading order in h. The classical turning points are 
Zl, Z2, and Z3. The WKB wave functions <j>i correspond to eqs. (3.7)-(3.12) 



boundary conditions that ~ (1 — z)^**^/^. The WKB wave function (3.12) has the near 
horizon expansion (pQ ~ ^6(1 ~ z)^*l'^l/^ + -66(1 ~ z;)*''^'/^. Thus when cj > we should take 
Bq = 0, and when w < we should take instead Aq = 0. The equation 

V6 = M5M4M3M2M1V1 (3.13) 

then provides a Bohr-Sommerfeld like quantization condition on uj. 

The matrices Mi are all well known. To go from a classically forbidden region to the 
classically allowed region, we use the standard WKB connection formula 

(3.14) 

Thus we find that Mi = M and M3 = M~^. To go from (p2 to cp^ or from cp^ to (p^, we make 
use of the fact that f{z)dz = — f{z)dz + J^^ f{z)dz: 





^2 I , M4 = ( ^ „ 1 , (3.15) 

where 

+ (3.16) 

The connection matrix M4 deserves closer scrutiny because the integral K is not always 
well defined. From eq. (2.21), we see that the potential term Vi (for the (p+ case) has the 
form 

Vi = -Qzzm = . (3.17) 

Thus Vi will have a simple pole where vanishes. From the form of Vq = Qzzfn^ — Z-^Zk, 
it is clear that Z^ will only vanish in a classically forbidden region where Vq > 0. Thus, the 
integral L will never be singular in this way. Let z = Oj be the locations of the simple poles 
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of Vi. Near aj, the integrand for K looks like ^{z — ai)^^. We regulate the singularity by 
taking a small semi-circular detour in the complex z plane. The detour introduces a factor of 
ibz7r/2, depending on the choice of detour above or below the singularity. This extra phase 
factor introduces a relative minus sign between the two nonzero entries of M4. We do not 
need to worry about the overall sign as it will not affect the quantization condition. 

We can be more precise about where has zeroes. From (2.19), is manifestly positive 
in the region < z < 1 if both w > and k < 0. Thus in this case, there will be no singularity 
to worry about. However, if u < there will in general be one such value z = ai and if w > 
and A; > 0, there will be either zero or two such values. 

Given the matrices Mi and the relation (3.13), we find the quantization condition on u 
looks in general like 

cosL(a;) + ice-2^(^)sinL(w) = 0. (3.18) 

In deriving (3.18), we made the implicit assumption that the turning points lie on the real 
axis. However, the quantization condition implies uj is not real, and if a; is not real, the turning 
points will in general not lie on the real axis. To get out of this apparent contradiction, we 
assume that |ce"2^| <C 1. Then the imaginary part of lo should be small, and we can use 
(3.18) to estimate it. We find the QNMs at ujn, n = 1, 2, 3, . . ., satisfy 



L(ReuJn) = T^in — 1/2) , Ima;„ 



dL/duj 



(3.19) 



In all our examples, it is also true that dL/duj > 0. Thus the sign of the imaginary part is 
determined by the sign of ce~^^. 

Specializing to the case for concreteness, there are several cases to consider. 

• When a; > and k < 0, the singularities at 2; = aj are absent, and we find for both the 
scalar and spinor the quantization condition (3.18) with c = —1/4. These QNMs lie in 
the lower half plane. 

• For the scalar when u; < 0, the change in the WKB boundary conditions at the horizon 
leads to the relation (3.18) with c = 1/4. These QNMs lie in the upper half plane. 

• For the spinor when a; < 0, there is a pole in Vi at z = ai. Deforming the contour to 
avoid the pole adds a phase factor TTi/2 to the integral K. Thus while c = 1/4 as in the 
scalar case, e~^^ < is negative. These QNMs lie in the lower half plane. 

• We may also consider the case where cj > and k > 0. In this case, Z/^ may have no 
zeroes or two zeroes in the region < 2; < 1. If there are no zeroes, we reduce to the 
cj > and k < case. If there are two zeroes, then e~^^ > and the QNMs are still 
in the lower half plane. 

At first order in h, the WKB approximation already works pretty well. Below are tables 
comparing the first order WKB against numerical integration of the Dirac equation. For 
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spin up 




spin down spin down 




Figure 10. Normal modes obtained by the Bohr- Sommerf eld quantization, marked by red dots. The 
parameters are m — 2, k = b, and /i^ — 25. Here \S\ — \aa\ is the source (denominator) of the Green's 
function. 



m 



2 and = 25 we have the following results when k = ±5, illustrating the spin splitting: 



/c = -5 


numeric 


WKB 


-7.90 - 0.00516i 
-4.23 - 0.00597i 
-1.42 - 0.00545i 
0.574 - 0.0390i 


-7.99 - 0.00504i 
-4.30 - 0.00587i 
-1.47- 0.0052i 
0.576 - 0.041i 



A; = 5 


numeric 


WKB 


-6.40 - 0.00398i 
-3.04 - 0.00476i 
-0.504 - 0.00471i 


-6.48 - 0.00383i 
-3.10 -0.00471i 
-0.540 - 0.0044i 



(3.20) 



We have presented the same information graphically in figure 10. 
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Figure 11. Quasibound states in tlie complex cj plane. The red dots are the quasinormal modes 
obtained by the generalized WKB formula. The parameters are m — 2, k — Q. Two values of fiq are 
considered. 



For m = 2, /c = 0, we also present results for /ig = 50 and = 75. 





Ptg = 


= 50 




numeric 


WKB 


-26.35 


-0.172i 


-26.77- 


- 0.135i 


-20.84 


- 0.248i 


-21.25 


-0.18i 


-16.29 


-0.281i 


-16.66 


- 0.19i 


-12.41 


- 0.292i 


-12.74 


- 0.18i 


-9.03 


- 0.290i 


-9.34 


-0.17i 


-6.10 


- 0.278i 


-6.37 


-Q.m 


-3.57 


- 0.255i 


-3.80 


-0.14i 


-1.44 


-0.218i 


-1.63 


-Q.l2i 







Aig = 


= 75 




numeric 


WKB 


-45 


98 


-0.214i 


-46.50 - 


- 0.165i 


-39 


14 


- 0.313i 


-39.64 - 


- 0.227i 


-33 


41 


-0.36H 


-33.87 


- 0.24i 


-28 


43 


- 0.385i 


-28.87 


- 0.24i 


-24 


02 


- 0.394i 


-24.42 


- 0.23i 


-20 


05 


- 0.394i 


-20.42 


- 0.22i 


-16 


45 


- 0.388i 


-16.80 


-0.21i 


-13 


18 


- 0.377i 


-13.51 - 


- 0.195i 


-10 


20 


- 0.362i 


-10.50 


-0.18i 


-7 


50 


- 0.343i 


-7.77 


-0.17i 


-5 


06 


-0.319i 


-5.31 


-0.15i 


-2 


90 


- 0.288i 


-3.11 


-0.13i 


-1 


05 


- 0.240i 


-1.23 


-O.lli 



(3.21) 



The results are presented graphically in figure 11. 

We will stop at first order in the WKB approximation, but there are some interesting 
complications that appear at second order. There exist second order poles in Vk^m at z = 0, 
1 and ai that require a more careful consideration. Regardless of the background, the second 
order poles near z = ai have the universal form 



3/4 

{z - ai) 



2 ■ 



(3.22) 
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In our particular case, at the conformal boundary, we find that 



Vk,m = '^\, '^^ + 0{z-'). (3.23) 

At the horizon, we have instead 

The issue with second order poles in the potential is that the naive WKB wave functions do 
not have the proper scaling behavior near such points. Langer [34] proposed a modification, 
justified by a rescaling of the wave function and redefinition of z, that boils down to adding 
by hand a term of the form \ {z — Zs)~'^ to the potential for each second order pole Zs- These 
Langer modification factors are second order in our h expansion, and thus we have neglected 
them. 



4 Discussion 



We have tried to show in this paper that spin affects in important ways the physics of holo- 
graphic constructions involving fermions. The pole of the field theory fermionic Green's 
function obeys a Rashba type dispersion relation, indicating the importance of spin orbit 
coupling in the gravitational side of the construction. The spin also plays an important role 
in determining the lifetimes of quasiparticles. Without the spin corrections in the WKB ap- 
proximation, the fermionic dispersion relation would have an imaginary part of the wrong 
sign! That spin plays an important role in the bulk is ironic given that in a 2+1 dimensional 
field theory dual, we argued in the introduction that the fermions should be both massless 
and (hence) spinless. 

This work leaves open several issues that we would like to return to at some point. 
One is a more thorough exploration of the parameter space of our model. We would like to 
understand better the qualitative difference between bulk fermions with large masses mL > 1 
and small masses mL < 1 illustrated by figures la and 4. While the gravity explanation is 
related to the absence of a rest frame for a relativistic particle, the field theory interpretation 
is less clear. The limit w — t- played an important role in previous works on the subject 
(see for example [9]), and we would like to see what our WKB formalism predicts for the 
magnitude of the imaginary part of the dispersion relation and the corresponding lifetime of 
the quasiparticles. (How does including h corrections change the results of [17]?) Third, it 
would be interesting to consider the T — )• limit in more detail. 

It would also be interesting to work with a wider variety of backgrounds, for example 
the electron star or the holographic superconductor. We believe the qualitative nature of 
our story involving spin-orbit coupling and quasinormal mode placement will not change, but 
there may be other interesting effects. For example, with the d-wave superconductor studied 
in [22], it was precisely this spin-orbit coupling which gave rise to Fermi arcs. 
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A Analytic results from Heun polynomials 

In the AdS/CFT correspondence, some perturbation equations can be written in terms of the 
Heun equation, which has four regular singularities. Under certain conditions, the solution 
of the Heun equation is a polynomial, which can help us to obtain an exact solution to the 
Green's function. The Heun differential equation is 



where a+/3+l = 7+5+e [35]. The regular singularities are at x = 0, 1, a, and oo. The solution 
to this equation is called the Heun function, if it is regular at both x = and x = 1 (assuming 
a > 1). We denote the Heun function by HeunG(o, Q; a, /3, 7, 5; x), which is symmetric in a 
and /3. The Heun function can be written as a power series HeunG(x) = Yl'^o'^r^'^ where 
the coefficients satisfy a three term recursion relation [35]: 

(r - 1 + a)(r - 1 + p)cr-i - {r[{r - 1 + + a) + aS + e] + Q}cr + (r + l)(r + 7)ac^+i = . 

Given this recursion relation, it is clear that if a or /3 = —n and if Cn+i = 0, then HeunG(x) 
is an nth order polynomial.^ In general, the condition c„+i = is an (n + l)th order algebraic 
equation in the parameters of the Heun function. Since we need to solve an (n + l)th-order 
algebraic equation to obtain an nth-order Heun polynomial, we usually cannot obtain explicit 
solutions when n > 4. (For comparison, the series expansion for the hypergeometric function 
2-Fi(a, /3, 7; x) is defined by a two term recursion relation, and the polynomial condition is 
just a = —n.) The set of solutions c^+i = may include unphysical regions of parameter 
space, for example regions where the charge or momentum is imaginary. We will need to 
further restrict the solution set. 



^There are other types of "Heun polynomials" [35]. For example, another solution to the Heun equation 
is 'HeunG(a, {aS + e)(l — 7) + Q; a + 1 — 7, /3 + 1 — 7, 2 — 7, 5; x), which indicates that the HeunG on 
the right-hand side is a polynomial of order n when a + 1 — 7 = —n. In this work, other polynomials give 
unphysical results or the same results as the case a — —n. 




(A.l) 



Setting 




(A.2) 
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the two coupled equations (2.8) for ui and U2 can be reduced to a single second-order ODE. 
To avoid the awkward square root, we define u± = ui±U2, and then the Dirac equations are 



i{uj + qAt) 
OzU+ ^i^ 

i{u + qAt) 
OzU- H u- 



m + ikz 



zVf 
m — ikz 



-U- 



-u+. 



(A.3) 
(A.4) 



The equation for n_ is given by^ 



u_ + 



+ 



m 



2/ z{m — ikz) 



u_ + 



(w + qAtf 



P 



2 

Z^f 



m 

T2 



i{uj + qAt 
f 



m 



f 

— 1 

2/ z{m — ikz] 



+ 



iqA^ 
f 



/ 
= 0. 



(A.5) 



We need to solve this equation with in-falling wave boundary condition at the horizon, and 
then plug n_ into eq. (A.4) to obtain n+. 

For the massless fermion in AdS^, we obtain a second-order ODE with four regular 
singularities as follows:^ 



f 

u_ + -u_ + 



{uj + qAt 
P 



fc2 if{uj + qAt) iqA[ 
f 2/2 + / 



U- 



0, 



(A.6) 



where f = 1 — z^ and At = 1 — z'^. The solution is 

n_ ~ (Z + l)-/4+l/2(^ _ ;l)-W4+1/2(^ ^ .)-g/2+l/2~.;/4(^ _ .)-3/2+9/2+a./4 

X HeunG 



1 3-{l-i)uj ik'^ 
2' 2 T 



3 — a; S + iwS — (1 — i){z + 1) 
9;2, g, , ■ 



2 z 



We can see that when 

w = 2n + 3 — 2g , ?i = 0, 1, 2, • • • , 
we can obtain an nth-order polynomial solution 

(l-i)(z + l) 



HeunG (x) = Pn = c^x'^ 



r=0 



2{z-i) 



(A.7) 
(A.8) 

(A.9) 



if k satisfies an (n + l)th-order algebraic equation. 
The zeroth-order polynomial is 



Po 



1 , when k'^ = 6 — 4q . 



(A.IO) 



^Note that this equation can be put in Schrodinger form and gives an alternate starting point for a WKB 
approximation. The corresponding Schrodinger potential is complex which makes phase integral methods 
substantially more involved. 

^While we have focused on AdS4 backgrounds in the bulk of the paper, with an appropriate choice of 
gamma matrices [29], eq. (2.8) also holds for AdS^. 
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With this exact solution, we can obtain the Green's function 

4-2q- V6-4g 



G{u = 3-2q,k= v^6 - Aq, 



(A.ll) 



4-2q + y/6-iq 

For the other sign of k, G{—k) = —1/G{k). We can see that Im(G) > is always satisfied (if 
k is real). The first-order polynomial is 



Pi 



1 



2lo -2iuj -6 + ik"^ + 4g 



iuj + 2, 
The Green's function 



when k^ = lb-Qq± yjAq^ - 20q + 33 . (A.12) 



G{uj = 5-2q,k = yi5-6q± W _ 20q + 33, q) 
can be exactly expressed. If we take the plus sign, for example, the result is 

4g+(7-2g)W-20g + 33 

4g + (7 - 2g)yV 



15 
15 



j3- 
(3- 



y/Aq^ - 20q + 33) J 15 - 6q + ^/Aq^ - 20q + 33 



(A.13) 



(A.14) 



20q + 33 + (3 + ^4g2_20(7 + 33)y^l5 - 6(7 + y^Aq^ - 20q + 33 

To obtain the second-order polynomial, we need to solve a third-order algebraic equation for k. 
There are higher-order polynomials, and the results are more complicated. The denominator 
of the above two Green's functions are non-zero for all real q. It is unlikely that we can solve 
for a Fermi momentum kp in this way. 

The exact results can be used to check the numerical program. For example, some exact 
results are 



G{oj 
G{uj 



l,k = 

l,k = 
-l,k 



V2,q = l) 

V6,q = 2) 
= 0,q = 3) 



(3- 
59 



2V2)i, 
- 24\/6 



(A.15) 

(A.16) 
(A.17) 



where the last one is consistent with G{uj,k = 0) = i for massless spinors [9]. As more 
examples, we list the numerical value of the exact result and the numerical result in some 
special cases as follows: 



w 


k 


Q 


exact 


numerical 


7 


\/21 + ^/57 


-1 


0.21336115i 


0.21336116i 


5 







0.16186619i 


0.16186620? 


3 




1 


0.10121094i 


0.10121094i 


1 




2 


0.04244923i 


0.04244924? 





23/4 


5/2 


0.04177353i 


0.04177353i 



There is a hidden algebraic structure in the Heun equation. A representation of the SU(2) 
algebra is 

S+ = z^d,-2sz, (A. 18) 

S- = -d,, (A. 19) 



So = zd. 



(A.20) 
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where [S+jS-] = 25z, [ScS'i] = ±S±, and S^ = s(s + 1). If we consider the fohowing 
Hamiltonian problem Hy{z) = Ey{z) with 



i,j=-,0,+ i=+,0, 
i>j 

we obtain a differential equation with polynomial coefficients [36]: 

P4 y"{z) + Ps y'{z) + P2 y{z) = Ey{z) , (A.22) 

where 

P4 = a^+z"^ + a+oz^ + (aoo - a+_)2;^ - a^-z + a — , (A. 23) 

P3 = 2(1 - 2s)a++z^ + [(1 - 3s)a+o + b+]z^ + [2s(a+_ - aoo) + ooo + bo]z (A.24) 
+ sao- - 6- , 

P2 = 2s(2s - l)a++2;^ + 2s(sa+o - 6+)^ + s^aoo - sbo ■ (A.25) 

There is a correspondence between the Heun equation and the spin system by the following 
identification: 

a-i-i- = , a+o = 1 ) CL = , (A. 26) 

aoo - a+- = -(1 + a) , ao- = -a , (A. 27) 

2s(a+_ - aoo) + aoo + &o = -(7(1 + a) + 5a + e), (A.28) 

sao- — b- = aj , s^ooo — sbo — E = Q , (A. 29) 

6+-3s = Q + /3, 2s{s-b+) = aP. (A.30) 

Solving eq. (A.30) gives s = —a/2 or s = —13/2. If the total spin s is an integer or half- 
integer, the H in eq. (A. 21) has a finite dimensional representation. Therefore, the eigenvalues 
satisfy an algebraic equation by diagonalizing H. This algebraic equation is equivalent to the 
condition for the existence of a Heun polynomial solution. If s is not an integer or half-integer, 
the Hilbert space is infinite dimensional. 

The interpretation of why there exist some exact solutions to the Green's function is 
as follows. When m = 0, in the three-dimensional parameter space {oj, k, q), there are 
discrete one-dimensional lines (subspaces) that are labeled by n = 1, 2, • • • . On these lines, 
the Hilbert space (accidentally) becomes finite dimensional. So the eigenvalues satisfy an 
nth-order algebraic equation. When n < 4, we can explicitly solve the algebraic equation. 
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